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1. Introduction 

In the present state of experimental particle physics, our strongest clues to the nature 
of physics beyond the standard model are the various fine tuning problems that are revealed 
by the standard model fit to experimental data. Surely the most pressing of these are 
the fine tuning of the cosmological constant, and the gauge hierarchy problem. Next in 
importance come the various small numbers associated with the fermion mass matrices, 
and among these, the strong CP problem is (at least numerically) the most striking. Several 
mechanisms, of varying degrees of plausibility, have been invented to account for the small 
value of the QCD vacuum angle that is required to explain the observed bounds on the 
neutron electric dipole moment. In particular, it is often suggested that if only a massless 
up quark were compatible with the spectrum of hadrons, it would also provide a neat and 
economical solution of the strong CP problem. 

Recently a class of models which employed discrete symmetries and some modest 
dynamical assumptions were proposed [PJ to explain the gross features of the quark and 
lepton mass matrices. The simplest of these models automatically incorporated a discrete 
symmetry which guaranteed that the up quark mass was zero to all orders in perturbation 
theory. The discrete symmetry enforced an accidental anomalous {7(1) symmetry on all 
renormalizable terms in the Lagrangian. This symmetry remains unbroken in perturbation 
theory despite the fact that the discrete symmetry is spontaneously broken. 

The real and imaginary parts of the determinant of the quark mass matrix are a priori 
independent parameters in the standard model. In the standard model, setting either or 
both of them to zero is unnatural. The argument that the theory is more symmetrical 
for Re m u = Im m u = 0, and that the fine tuning is therefore natural, is spurious. The 
axial i7(l) "symmetry" of the standard model with one massless quark is anomalous, and 
is really no symmetry at all. In the absence of further constraints on very high energy 
physics we should expect all relevant and marginally relevant operators that are forbidden 
only by this symmetry to appear in the standard model Lagrangian with coefficients of 
order one. 

In some of the models of 0, it is natural to set the real and imaginary parts of the 
up quark mass matrix to zero at high energy. The accidental anomalous {7(1) symmetry 
of these models guarantees that both the real and imaginary parts are generated only by 
nonperturbative QCD processes. In section II of the paper we review the details of the 
mechanism which leads to the accidental anomalous {7(1) symmetry, and argue that this 
is essentially the unique way to generate such an accidental symmetry. 
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To evaluate the viability of these models, we were led to reexamine the controversy 
which has surrounded the question of whether a massless up quark is consistent with low- 
energy hadron phenomenology In section III of this paper we review the literature on 
this subject. Fitting low energy data to a first order chiral Lagrangian, one can extract 
the "low-energy quark masses," These should be distinguished from the quark mass 
parameters, mi of the QCD Lagrangian at high scale (say 1 TeV) || |§. The distinction 
between them is second order in the m's. In particular, \i u receives an additive contribution 
of order ™ d ™* where A. x $b ~ 1 GeV is some characteristic QCD energy scale, perhaps the 
fundamental scale 47rf n . Even if m u = 0, the parameter jj, u can be nonzero. Its value is 
\i u = (3 ™ d ™° where the dimensionless coefficient, (3, is estimated to be a number of order 
one. Unfortunately, with present technology we are unable to calculate it reliably starting 
from QCD. The low-energy data is compatible with m u = provided (3 ~ 2 which is of 
order one. Therefore, m u = seems to be a viable possibility. Kaplan and Manohar [ffl 
tried to avoid a calculation of the parameter (3 by extending the low-energy analysis to 
second order in the m's. This led them to find an ambiguity in the parametrization of the 
second order Lagrangian which prevented them from determining the value of m u . Several 
authors § § @ § § tried to resolve this ambiguity and to determine m u by adding 
more physical input. Our understanding is that this ambiguity cannot be resolved solely 
on the basis of low energy data. Therefore, we conclude that the possibility that m u = 
is still open. 

In section IV we take up the question of whether a vanishing up quark mass really 
solves the strong CP problem. In models with an accidental anomalous U(l) symmetry, 
both real and imaginary parts of the up quark mass are generated by nonperturbative 
QCD dynamics. The message of section III was that the real part so generated is, within 
our ability to calculate, compatible with low energy data. Is the same true of the imaginary 
part? This is a particular case of a general question first addressed in 1979 by Ellis and 
Gaillard flO ], and since studied by a number of authors [jll|] ||12|1 : if the QCD vacuum angle 
is set to zero at some large scale Aq, what will its low-energy value be? We argue that in 
the present context there is a very general operator analysis of this question. The analysis 
shows that CP violation in the high-energy theory can induce a low-energy Qqcd only via 
the action of a CP violating irrelevant operator O. In this case, the same operator will 
make a direct perturbative contribution to the neutron electric dipole moment of order 
AxSB Oqcd where a u is the low energy up quark mass described above and A y sb again 
denotes a typical QCD scale. In other words, if both the real and the imaginary part 
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of m u vanish at high-energy, nonperturbative strong CP violation will be smaller than 
perturbative contributions to the neutron electric dipole moment. We conclude that the 
models of |I| provide a phenomenologically viable and relatively economical solution to the 
strong CP problem, in a framework which may resolve many of the other puzzles of the 
fermion mass matrix in the standard model. 

In an appendix we present a warmup exercise for the calculation of direct contribu- 
tions to the neutron electric dipole moment in models with an accidental anomalous U(l) 
symmetry. 



2. m u = 0, Naturally 

CP violation in the QCD Lagrangian arises from the terms: 

£cp = -f-^ G^G^ - iq(Im M q )~f 5 q (2.1) 

where M q is the quark mass matrix. There is only one independent CP violating parameter 
in (§, 

9 = 9 + argdetM. (2.2) 

Strong interactions, through their ^-dependence, lead to an electric dipole moment for the 
neutron. An estimate of this contribution, using current algebra, gives [|13| 

V n = +3.6 x 10" 16 9 e cm. (2.3) 

The experimental bound on V n [ jT3"| , 

V n < 1.2 x 10" 25 e cm (2.4) 

implies 

9 < 10" 9 . (2.5) 

In pure QCD the fine tuning required to satisfy ( |2.5|) is natural. A new symmetry, CP, 
is obtained for 9 = 0. However, within the Standard Model, this fine tuning is not natural, 
since the standard model lagrangian with nonzero CKM phase S is not CP invariant even 
for 9 = 0. The experimental observation of CP violation in the neutral K system requires 
that S be nonzero. 
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The above parametrization gives the erroneous impression that one fine tuning, m u = 
0, sets the up quark mass to zero and solves the strong CP problem simultaneously. In 
fact, if we choose to absorb the vacuum angle into the quark mass matrix, we can instead 
parametrize QCD by the real positive values of Nf — 1 quark masses, and the real and 
imaginary parts of the up quark mass!. The bound on 6 can be reexpressed as 

Imm„ < 5 x 10" 12 GeV, (2.6) 

which is about fourteen orders of magnitude below the electroweak breaking scale. Given 
this bound, the imaginary part of the up quark mass makes a completely negligible contri- 
bution to hadron masses, and chiral lagrangian fits to pseudoGoldstone boson masses are 
in fact constraints only on the independent, CP conserving, real part of the up quark mass. 
The latter is only constrained to be about five orders of magnitude below the electroweak 
scale, and we will see that the data are consistent with it being equal to zero. 

In order to render the vanishing of both the real and the imaginary parts of m u natural 
(with all other quarks massive), a continuous U(l) symmetry 

u -> e ia u (2.7) 

must be imposed. It guarantees that the left handed field u couples only to the gauge fields. 
Like the Peccei Quinn symmetries of axion models, this U(l) symmetry is anomalous and 
hence not an exact symmetry of the full theory. Therefore, it is unnatural to impose it 
on a Lagrangian at a given energy scale, unless we can provide a reason why symmetry 
violating physics at much larger energy scales does not violate all conclusions based on 
the anomalous U(l). We conclude that, the value of neither the real nor the imaginary 
part of the up quark mass is natural in the standard model, although the imaginary part 
is certainly the worst offender. 

We should remark that in SUSY theories, the non-renormalization theorems allow us 
to set the coupling of the up quark to the Higgs to zero. In such theories a massless up 
quark might appear to be technically natural. (Clearly, the symmetry ( |2.7| ) acts then also 
on the u squark.) However, unless we arbitrarily impose the U(l) symmetry on the soft 
SUSY breaking terms (in particular a coupling of the u squark to a squark and a Higgs 

4 In order for the eigenvalues of Re M to have the conventional meaning of quark masses, it is 
necessary to have Im M proportional to the unit matrix. The difference between Re M and the 
quark masses in our convention is 0{6 2 ) and therefore negligible. 
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breaks the symmetry) radiative corrections would lead to couplings of u to the Higgs and 
to an up quark mass. Since we do not expect the full theory to respect the anomalous 
U(l), the non-renormalization theorems do not help. 

The only natural way to ensure an anomalous U(l) symmetry is to make it an acci- 
dental symmetry. That is, it follows from another symmetry plus renormalizability. All 
terms which violate it are then irrelevant in the renormalization group sense, and will 
have negligible effect if the physics responsible for them is pushed off to a sufficiently high 
energy scale. Thus our strategy is as follows: We impose an anomaly free symmetry H, 
which can be continuous or discrete, and require that the most general renormalizable H 
invariant Lagrangian exhibits the accidental U(l) symmetry Q2.7| ). 

It is clear that in order to achieve this state of affairs, the symmetry H must act 
differently on u and on c and t. (It can also act on the other light fields.) It therefore 
satisfies the definition of a horizontal symmetry. As such we can immediately apply some 
of the general analysis of horizontal symmetries presented in . 

If H acts only on «, it is anomalous. To cancel this anomaly other fields in the 
theory must be in a complex representation of H. Since they are massive, H must be 
spontaneously broken. 

We now examine how H can be spontaneously broken. The expectation value of the 
single Higgs in the standard model cannot break H [Q. More precisely, a subgroup of 
U{l)y x H isomorphic to H is always unbroken (U(1)y is hypercharge) . Similarly, in 
a two Higgs version of the standard model with natural flavor conservation (such as the 
supersymmetric standard model) H cannot be broken by the expectation values of the two 
Higgs fields <p u and <pd [0- To show that, note that in the absence of a 4> u 4>d term in the 
superpotential and in the soft breaking terms, the renormalizable theory is invariant under 
the anomaly free U(l) symmetry 

u -> e~ Sia u 
e ia d 

s^e ia s (2.8) 
b -> e ia b 
<Pd - e~ ia $d 

(all other fields are invariant) which guarantees the accidental U(l) (|2.7|) . However, its 
breaking by (4>d) leads to an unacceptable Goldstone boson. With the 4> u (j>d term in the 
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Lagrangian we can again use U(1)y to make both Higgs fields invariant. Therefore, we 
must add more fields to the theory. 

The simplest way to organize the analysis is to integrate out the extra fields, which 
we assume to be heavier than the weak scale, and to consider non-renormalizable terms 
added to the standard model Lagrangian. Typical non-renormalizable terms which could 
avoid the previous no-go theorem are 

Qtf^)\Q^*(*g)' (2.9) 

where M is the scale of the fields which have been integrated out. Such terms explicitly 
break the anomaly free continuous symmetry ( |2.8| ) to a discrete subgroup in the absence 
of a 4>u<fid term in the Lagrangian. We can also add a field, S, which is invariant under 
the standard model gauge group, and whose vacuum expectation value breaks H. Then, 
terms like n n 

Q<M (lf) + Q ^*(lf) (2 ' 10) 
can lead to the entries in the mass matrix. Rather low values of (S) and M are consistent 
with the constraints on flavor changing processes [|l|. Terms of the form ( |2.9| ) ( pTT0|) can 
be generated by integrating out massive fermions of masses of order M [[15| and lead to 
interesting mass matrices for the quarks. 

We conclude that the framework of |l| provides the unique natural way of enforcing 
the accidental U(l) symmetry, and hence m u = 0, to all orders in perturbation theory. We 
must now examine the question of whether a theory with such an accidental anomalous 
symmetry is phenomenologically viable. 



3. Microscopic and Macroscopic Mass Matrices 

In this section we will review the controversy that has revolved around the "experi- 
mentally determined" value of the up quark mass. Early estimates [[TB| of the pion and 
kaon masses in terms of quark masses in lowest order X-PT led to the conclusion that the 
up quark mass could not be much smaller than half the down quark mass. 

In quantum field theory, the parameters in an effective Lagrangian are scale dependent. 
A microscopic theory more fundamental than the standard model fixes them at some high- 
energy, e.g. of order 1 TeV. In order to relate these to low-energy parameters we must 
understand how the parameters renormalize in QCD. 
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To lowest order in chiral perturbation theory, the renormalization of quark masses 
is multiplicative and flavor independent. This is what enabled Weinberg [|16| to obtain 
renormalization invariant predictions for high-energy quark mass ratios in terms of low- 
energy meson masses. Once we go beyond lowest order xPT however, this is no longer 
true. In particular, the up quark mass receives an additive renormalization proportional 
to m* d m* via the mechanism of || 1| . Their result may be phrased as the statement that 
a renormalization group transformation from a scale A to A — AA (for A ^> ^ x sb ~ 1 GeV 
where the theory is weakly coupled) leads to 



where the mass matrix M q has complex eigenvalues (the vacuum angle 9 has been 
rotated into M q ). The second term arises from small instantons and the dimensionless 
constant (3 depends on the number of colors N and flavors Nf. 

This term also depends on the regularization scheme through our choice of integrating 



It is of order exp(— pjxy) <C 1. Therefore, even though rrii(X) is ambiguous, the limit 
lim^oo mi(X) is well defined and M q at high energies is not ambiguous. 

In renormalizing between two high-energy scales (e.g. 1 TeV and 100 GeV) the additive 
renormalization is accurately calculated in the dilute instanton gas approximation, and 
is very small because it is nonperturbative in the small high-energy value of the QCD 
coupling. When the QCD coupling becomes strong this is no longer the case and the 
best we can do is to estimate the additive contribution to the low energy quark mass as 
Therefore, if m u = 0, fi u ~ T^fj • ^his estimate is compatible with the low- 
energy quark mass ratios extracted from lowest order chiral perturbation theory. A more 
detailed comparison requires a nonperturbative solution of QCD, and might in principle 
(though probably not in practice since the effect pointed out in || and @ is invisible in 
the quenched approximation) be extracted from a lattice calculation. 

To conclude this section, we emphasize that the high energy values of the quark masses 
are, for all practical purposes, unambiguously determined by potentially observable QCD 
Green's functions. For example, in the QCD sum rule approach to hadron phenomenology, 




(3-1) 



over the instanton scale size p from A 



to (A — AA) . This ambiguity is very small. 
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it is precisely these high energy values which are related to low energy hadron parameters!. 
Unfortunately, it is difficult to estimate the uncontrollable systematic errors in sum rule 
calculations. As far as we have been able to determine, all such calculations are compatible 
with a vanishing value for the high energy up quark mass. 

3.1. First order analysis 

To make these considerations more precise, let us outline the procedure for actually 
extracting the high-energy quark masses from low-energy data. The first step is to fit 
low-energy data to the low-energy chiral Lagrangian (we use the notations of reference [JL7[ 

F 2 

C= -^-Trid^Ud^ + 2B(Mrf +M^U)}. (3.2) 

Using SU(3)l x SU(3)r, A4. which is in the (3,3) representation, can be brought to the 
form 

M = I // d I (3.3) 

where all the parameters are real. The mass parameters \ii can be interpreted as the 
low-energy values of the high-energy quark masses mj. The result of the fit is 

0<1O- 9 , ^-0.57, ^-0.05. (3.4) 

fJ-d Us 

Now let us expand the ^'s in the short distance m's. Using the SU(3)l x SU(3)r 
transformation laws of the m's we find that for three flavors 

777 * 771 * 

A XSB 
777 * 777 * 

\x d = Pim d + /3 2 -r^ + °( mS ) (3-5) 

A XSB 
A XSB 

where (3\ and fa are dimensionless coefficients. The /x's in ( |3.5| ) are complex and can 
be brought to the form ( |3.3|) using an SU(3)l x SU(3)r transformation. The additive 



5 This implies that if our assertions about the inequality of low energy and high energy quark 
mass ratios are correct, then values for the quark masses which are extracted by combining sum 
rules with chiral lagrangian results for mass ratios are incorrect in principle. 
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renormalization proportional to j3 2 is the strong coupling version of the instanton term in 
( ft.lp and it represents the breaking of the axial U(l) by the anomaly. In terms of the 
matrix A4 and the underlying quark mass matrix M q , equation (|3.5|) can be written as 

M = PiM q + f3 2T ^—detMl—^M q + £>(M 3 ). (3.6) 
A x sb q M q M\ ql 

The dimensionless coefficients (3\ and j3 2 cannot be found without a strong coupling 
calculation in QCD. On general grounds we expect them to be of order one. They both 
suffer from ambiguities resulting from the regularization scheme. The only sense in which 
any of them is small is in the large N approximation, where (3 2 is of order 



Assuming ^ <C 1, (|3.5|) leads for real m's to 



m d \i d Pt A xSB \ A xSB J 

Therefore, — cannot be determined without knowledge of |j§. If the dimensionless ratio 
■fl is near 2, the up quark can be massless. Even if this ratio is not near 2, the ambiguity 
in — is significant. 

3.2. Second order analysis 

The previous discussion can be criticized on the basis that it is a first order analysis 
in the //'s but it includes some second order contribution in the m's. To make the analysis 
consistent one should include all second order terms in the m's and hence, go to second 
order in the /i's. The most general potential of second order in M q can be written, using 
X = 2BM q , as 

V(U) = j^Trxtft + ri[Trx t ^X t f/- (Tr^t/) 2 ] j + h.c. 

+r 2 [Tr X i Ux t U + (Trx^) 2 ] + h.c. ( 3 - 8 ) 
+r 3 (Trx t f/)(Trxf/ t ) + 0(M 3 ). 

(The dimensionless coefficients rx, r 2 and r 3 are linear combinations of Lq, L-j and L§ of 
|j~7l1.) Remembering that M q is in the (3, 3) representation of the SU(3)l x SU(3)r flavor 



symmetry, we learn that the terms multiplying t\, r 2 and r 3 are in the representations (3, 3), 
(6, 6) and (8, 8), respectively. Since the two terms in the curly brackets in (|3.8| ) transform 



9 



as one irreducible representation, there must be a redundancy in the parametrization. 
Indeed, using the identity of 3 x 3 matrices 

det A - -A[(TtA) 2 - TrA 2 } + A 2 TtA - A 3 = (3.9) 

applied to A = x^U we find 

det ^Tr^ptft + ^TrxWU - ^(Tr X + C7) 2 = 0. (3.10) 

Therefore, the shifts 

Sa . . + . 1 
X^X + -^detx T — 

^ X T (3.11) 

f*i — »ri + a 

with an arbitrary complex number a change the potential ( ^.8| ) by terms of order Mjj, 
which are neglected. 

Because of this redundancy in the parametrization of the effective Lagrangian, the 
mass matrix M q cannot be determined by fitting the experimental data to (|3.8| ). At best, 
we can determine M q up to an arbitrary shift by ^§(det Mj)-^. 

The ambiguity ( |3.11| ) is reminiscent of the expression for M. in terms of M q ( |3.6| ). 
The functional similarity stems from the fact that the axial U(l) is broken preserving only 
the SU(3)l x SU(3)r symmetry. However, these two effects are different. Equation (|3.6Q 
relates the short distance mass M q to the long distance mass M. defined by a fit to a first 
order Lagrangian. The non-linear term in this equation expresses a renormalization effect. 



The ambiguity (|3.11|) expresses a redundancy of the second order low-energy description 
which prevents us from using only the low-energy chiral Lagrangian to determine M q . 

The distinction between these two effects is more clear in the "trivial" two flavor case. 
The identity analogous to ( |3.10| ) in this case is 

TrV' t C/ = detV t Tr- T C/ t (3.12) 

where ip is an arbitrary 2x2 matrix. It shows that the matrix x m the potential Trx"^ U+c.c 
is ambiguous. First, it is clear that only the eigenvalues of x are physical. Second, using 
equation ( |3.12| ) the transformation 

x ^ x + ^_ det ^ti_ ( 3 . 13 ) 
10 



leaves the potential invariant and allows us to set one of the eigenvalues of x to zero!. The 
non-perturbative contribution to M. arises at first order in M q . It is the second term in 

M = PiM q + (3 2 det M\ — + 0(M O 2 ) . (3.14) 

q M q M\ q 

Clearly, the redundancy in the parametrization has nothing to do with the change in the 
strength of this effect (changing fa)- The only common thing about the ambiguity ( |3.13j ) 
and the additive renormalization ( |3.14|) is that they both allow for m u = to be compatible 
with low energy data. 

Some authors || || [0 have suggested the use of more input about low energy hadron 
physics to resolve this ambiguity and determine M q . We would like to make clear what 
bothers us about these attempts. These authors fix the ambiguity by insisting that the 
coefficients x m Tr\U^ in the chiral Lagrangian are proportional to the bare quark masses, 
i.e. by setting what we have called \i u oc m u . They then attempt to compute the other 
coefficients in the Lagrangian (implicitly) with this choice of Kaplan-Manohar "gauge." 
However, their computations make no reference to this particular choice. We could equally 
well consider them to be determinations of the couplings in some other "gauge." And since 
their computations involve relatively low energy hadronic physics there is no reason to think 
that they have anything to do with the special parametrization in which the low energy 
and high energy up quark masses are proportional to each other. Other approaches of 
resolving the ambiguity use the large iV expansion. While this may be a valid approach, 
one must beware of considerations that are justified only in this limit. The whole issue is 
of higher order in jt and therefore cannot be settled at the leading order in this expansion. 

Finally, we note that the argument is sometimes made that the consistency of the 
quark masses extracted from first order chiral perturbation theory analysis of baryons 
and mesons, suggests that the scale dependent additive renormalization which we have 
discussed, is small. That is, the additive renormalization should give an effective upquark 
mass proportional to mdfn s in both the baryon and meson effective lagrangians, but the 
coefficients would not be the same. We have examined the extraction of quark masses 
from baryon masses Jl8j || and believe that it suffers from large uncertainties. The results 

3 /2 

are sensitive to nonanalytic corrections (oc m q ), and to model dependent calculations of 



6 Note that therefore, in the two flavor case there cannot be any CP violation in the first order 
chiral Lagrangian. 
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electromagnetic contributions. They are not precise enough to rule out the possibility that 
m u = 0. 

We conclude that at the level of precision (order of magnitude) of nonperturbative QCD 
calculations available to us at present, low-energy phenomenology is completely compatible 
with a vanishing value of the high- energy up quark mass. 

Only a nonperturbative calculation in QCD can prove or disprove the phenomenolog- 
ical viability of m u = 0. Therefore, in view of the recent progress in numerical methods 
in lattice gauge theory, we would like to encourage a detailed analysis of the possibil- 
ity of a massless up quark by these methods. We emphasize however that the additive 
renormalization (|3.1| ), vanishes in the quenched approximation. 

4. Does m u = Solve the Strong CP Problem? 

As we have seen in section II, models of the type studied by [0 can incorporate discrete 
symmetries which guarantee that both the real and imaginary parts of the up quark mass 
are zero (at 1 TeV) while CP violation in the Cabibbo-Kobayashi-Maskawa matrix is 
allowed. The discrete symmetries enforce an accidental anomalous U(l) symmetry on the 
Lagrangian. If all physics not explicitly included in the Lagrangian comes from energy 
scales substantially higher than 1 TeV then this accidental symmetry guarantees that the 
only nonvanishing contributions to the up quark mass come from nonperturbative QCD 
effects. 

We have argued above that in our present state of impotence with regard to precise 
nonperturbative calculations in QCD, the real part of the up quark mass generated in 
this manner appears to be consistent with hadron phenomenology. The purpose of the 
following discussion is to determine whether an analogous conclusion can be made for the 
imaginary part of the up quark mass. 

This task is made easier by the observation that violations of the anomalous U(l) 
symmetry are only significant at scales below 1 GeV. Thus if we construct an effective 
field theory for scales just above 1 GeV, it must obey the U(l) symmetry. As a conse- 
quence, there will be no CP violating terms in the renormalizable part of this effective 
Lagrangian. The imaginary part of the low-energy up quark mass will thus be generated 
by a combination of nonperturbative QCD effects and CP violating irrelevant operators 
in the effective field theory. Let Od be such an operator, of dimension d. It will give a 
contribution to the imaginary part of the up quark mass of order Cf^(^f^) d " 4 , where 
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M d-4 is the coefficient of Od in the effective Lagrangian, M is the heavy scafe at which 
CP violating effects originate (e.g. the mass of the W boson), and A x sb is, as always, a 
typical QCD scale of order 1 GeV. The consequent contribution to the neutron electric 
dipole moment is of order C^^(Afp) d - 4 . 

XSB 

The operator Od also makes a direct contribution to the neutron electric dipole mo- 
ment of order C-^ — ( A ^| B ) d ~ 4 , which is larger than that coming from the "induced strong 

* A 2 

CP violation" by a factor of order xSB . Thus, if the direct contributions to the neutron 

J morris > J 

electric dipole moment coming from perturbatively induced CP violating irrelevant opera- 
tors are within experimental bounds, there will be no strong CP problem. Perturbative CP 
violating effects might put significant constraints on particular models with new physics 
at the TeV scale. We will not study these direct contributions in detail here (except for a 
short discussion of the standard model in the appendix). What we have shown is that if 
these constraints are satisfied, we need not worry about strong CP violation. 
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Appendix A. Direct Contributions to the Neutron Electric Dipole Moment 

Here we present an operator analysis of direct perturbative contributions to the neu- 
tron electric dipole moment in the standard model with an accidental U(l) symmetry. 
This is a warmup for a full calculation in models of the sort studied in . 

The lowest dimension CP violating operators that might be relevant here are the 
dimension five chromoelectric dipole moments of the light quarks. In models with an 
accidental U(l) symmetry, no dipole moments involving an up quark are allowed in the 
effective Lagrangian above the QCD scale. It is also worth reiterating |TI| that all dipole 



moments are suppressed by factors of light quark masses because they violate the non- 
abelian chiral symmetries of massless QCD. As a consequence their contributions to the 
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neutron electric dipole moment are a priori smaller than chirally allowed dimension six 
operators, including the three gluon operator PH , and the axial polyp [T^] operators: 

There is a similar term, with a photon field i* 1 ^ replacing the gluon field G^ v in ( |A.1|) . 
In the standard model, the largest contributions to quark electric dipole moments at a 
given scale actually come from "fusing" together one of these polyp operators and a quark 
mass term via the exchange of a gluon right at the infrared cutoff scale. In many previ- 
ous analyses of CP violation, the electric dipole moment operators were considered to be 
more important than, or on equal footing with, the polyp operators. We believe that such 
arguments rest on the dubious procedure of sticking a "constituent" quark mass into an 
effective lagrangian. The purpose of an effective lagrangian calculation is to make a clean 
separation between physics at different scales, and in particular between weakly coupled 
high energy physics, and strongly interacting QCD. We do not believe that it is consistent 
with the "rules of the game" to calculate the coefficients of an effective lagrangian (with- 
out the use of a computer) at scales at which QCD is strongly coupled. Thus, constituent 
quark masses should not appear as coefficients in an honest effective lagrangian. Rather, 
we should work with operators normalized at a scale where QCD is still weakly coupled and 
estimate their hadronic matrix elements by dimensional analysis or some more sophisti- 
cated hadronic "model." In the framework of a such a philosophy, the direct contributions 
of CP violating polyps are more important than chromoelectric dipole moments of light 
quarks. 

At the order of magnitude level then, the calculation of the neutron electric dipole 
moment in the models of [l|] probably reduces to the calculation of A3 for the various 
quarks. As a warmup exercise for this calculation we will present an estimate for the A3 
coefficient within the Standard Model. 

As a single VF-propagator cannot introduce CP violation, we are led to consider the 
diagrams of order g s o? w shown in Fig. 1-3. The dependence of such diagrams on the 
quark sector parameters is of the form 

3 3 



J2 E A(m)f2(m j )f 3 (m k )Im[V ui VlV jk Vl 

i,k=l j = l 



(A.2) 

33 3 v ' 

= E ^h( m i)fc{ m j)h{ m k) J E 6 ^rn^ikn- 
i,k=l j = l m,n=l 
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J is the CP violating invariant measure of Jarlskog [21]. As only left handed fields have 
charged current interactions, there can be no odd number of helicity flips in the loops, 
so the fi functions are functions of quark squared-masses, m 2 q . Assume that we had 
replaced f2(jnj) by a function independent of the quark masses. Since ^ ■ eij m = 0, such 
a quantity would not contribute to the coefficient of CP violating operators. Similarly, if 
either of fi(rrii) and fs(mk) were constant, there would be no contribution. Therefore, 
we can rewrite the above quantity by adding and subtracting contributions which do not 
contribute to CP violating terms. The result is: 



J{[fi(m 2 s ) - h{ml)][f 2 {ml) - / 2 (m t 2 )][/ 3 (mi) - / 3 (mg)] - (m* <-> ml)}. 



(A.3) 




Fig. 1: First Feynman diagram for induced 8 




Fig. 2: Second Feynman diagram for induced 8 
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(c) 



Fig. 3: Third Feynman diagram for induced 9 



Next, let us concentrate on the inner loop of diagrams 1(a) and 1(b). We assume that 
the momentum p entering this loop is small on the scale of mw (This assumption will be 
justified below.) Then such a loop contributes 



£1 



d 4 



i> 



(/ + p) 2 + m 2 w \l 2 + m 2 c l 2 + m 2 



= {m\ - ml) J d L 



((l + p) 2 + m 2 w ) {I 2 + ml) {I 2 + m 2 ) ' 



We next expand in powers of p. The zeroth order term does not contribute because of 
Lorentz invariance. The first order term is 



Li ~ (m 2 - m\)i> 



dH V 



(I 2 + m 2 w ) 2 {l 2 + ml)(l 2 + m 2 t ) ' 



(A.4) 



We can scale the integral by m 2 , assumed to be of order m\y. The integral is finite, 
and consequently L\ ~ j>. If mt had been smaller than mw, there would have been an 



additional suppression factor of order -^f 



Thus we can replace the inner loop with an insertion of a ^-operator, where p is the 
momentum entering this loop. In this approximation, the external loop of the diagram 
with the gluon on the left (Fig. 1) gives 



x 



i 



i 



" /m ((l-k) 2 + m 2 w ) 



1 



1 



(l 2 + m 2 s ) (l 2 + m 2 ) 
1 



((J + q) 2 + m 2 )(l 2 + m 2 ) ((I + q) 2 + m 2 )(l 2 + m\) 
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- (m s <-> m d ) 



(A.5) 



where q is the incoming gluon momentum and k is the incoming w-quark momentum. The 
diagram with the gluon on the right (Fig. 2) gives a similar contribution in which the 
coefficient of ^1^2^^ is modified by simultaneously exchanging rrid with m s and I with 
I + q. If we now add these diagrams and perform the explicit subtraction of terms with 
rrid <->• m s , we obtain: 

L 2 ~(m s -m d )(m b - m a )(m b - m d ) ' 



((I - /c) 2 + m^)((Z + g) 2 + ml) (I 2 + m 2 ) 



1 

x ■ 



((/ + q) 2 + m 2 )(l 2 + m 2 )((l + q) 2 + m 2 )(l 2 + m 2 ) ' 

(A.6) 

We are interested in the terms linear in both k and q. Thus we expand to first order in k and 
q. As Li is a highly convergent integral even when we neglect the momentum dependence of 
the W propagator, we were justified in assuming that the momentum entering the internal 
loop was small. We obtain a factor m^, multiplying a finite integral. Scaling the loop 
momentum by mj and expanding in (^-) 2 and (^-) 2 we get 

T < f dHl-kl-q^J 

2 ~<J (/ 2 + l) 2 (/ 2 + (^) 2 ) 2 (/ 2 + (^) 2 ) 2 ' 1 ' } 

This leads to@ 

A^(^) 2 J^L (A.8) 

Using this result, it is easy to estimate the induced 9. We can "fuse" the polyp 
operator with the up quark mass (which in our case is induced dynamically) to find 



A 2 2 

Im m u ~ (Re m u )A 2 ss A^ ~ J(^) 2 (Re m u ) x ^ s (A.9) 



which agrees with the estimate of \W] . Note that with the extra suppression factor coming 

2 

from L\ ~ ~ 7 i~ r h when mt < rnw, our estimate (IA.9I) would agree with that of Shabalin 

m w 

ph. 



For the diagram in Fig. 3, the inner loop can be replaced with a gluon vertex (? S 7 M . It 
contributes to A3 with the same order of magnitude as diagrams (a) and (b). 
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